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Abstract 

The spectrum of two-particle bound states is investigated within a relativistic 
quantum-field model of interacting quarks and gluons confined analytically. The 
hadronization process of mesons and glueballs is described by using the Bethe- 
! Salpeter equation. Provided by a minimal set of physical parameters (the quark 

' masses, the coupling constant and the confinement scale), the model satisfactorily 

describes the meson ground-states, orbital and radial excitations in a wide range 
from 140MeV to 9.5GeV. The estimated values for the coupling constant and the 
lowest-state glueball mass are in reasonable agreement with experimental data. 
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^ 1 Introduction 

Conventionally, the observed colorless hadrons are considered as bound states of quarks 
and gluons under the color confinement of the QCD. The quark- ant iquark pair fits 
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' in a pattern, where they have quantum numbers of the most of mesons. The color 



O 
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confinement is achieved by taking into account nonperturbative and nonlinear gluon 



\Q • interaction and the coupling becomes stronger in the hadron distance [1]. The correct 



summation of the higher-order contributions becomes a problem. Besides, the structure 



i-£h ' of the QCD vacuum is not well established yet. 

The conventional QCD encounters a difficulty by defining the explicit quark and 
gluon propagator at the confinement scale. Particularly, the Schwinger-Dyson equation 
^ ■ is used to obtain an effective quark propagator but this implies additional assumptions 
on the behaviour of the gluon propagator and quark-gluon vertices and results in elab- 
orated numerical evaluations [2]. 

A number of phenomenological approaches is devoted to the hadron spectroscopy. 
Some of them are adapted to the specific sector of heavy hadrons and introduce many 
parameters, including nonphysical ones. The 'potential' models use a stationary nonrei- 
ativistic Schrodinger equation with an increasing potential to describe the hadroniza- 
tion process of few-body systems. However, the binding energy in hadrons is not 
negligible, e.g., the observable two-quark bound states have masses much heavier than 
the predicted combined masses of the 'current-quarks' and it requires a relativistic 
consideration. 

It seems reasonable to develop a simple relativistic quantum field model of interact- 
ing quarks and gluons and consider the formation and spectra of the hadrons by using 
the Bethe-Salpeter equation. There exists a conception of the analytic confinement 
based on the assumption that the QCD vacuum is realized by the vacuum self dual 
gluon field which is the classical solution of the Yang- Mills equation [3, 4]. Accord- 
ingly, each isolated quark and gluon is confined and the corresponding propagators are 
entire analytic functions in the background gluon field. However, direct use of these 
propagators leads to long and complicated estimations. 
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In our previous investigation, we have clarified the role of the analytic confinement 
in formation of the hadron bound states by using simple scalar-field models [5]. We 
have considered "mesons" composed from two "scalar quarks" glued by "scalar gluons" 
by omitting the spin, color and flavor degrees of freedom. Despite the simplicity, these 
models resulted in a quite reasonable sight to the underlying physical principles of 
the hadronization mechanism and could give a qualitative description of the "scalar" 
mesons [6]. 

Below we extend the consideration by taking into account the spin, color and fla- 
vor degrees of freedom for constituent quarks and gluons within a simple relativistic 
quantum field model. This allows us to obtain extended and accurate results not only 
for the meson ground states but also for the orbital and radial excitations as well as 
the coupling constant and the glueball lowest state. 



2 The Model 

Consider a Yukawa-type interaction of the quarks ^^(x) and gluons 4> a ^(x) confined 
analytically. The model Lagrangian reads [7] 

c= (k[s-%%)+1 {<t>;[D-xJt)+g {K[^r%%<f>;)+g r bc , (i) 

where g - the coupling constant, F^ v = <9 M 0£ — <9„0£ and {a, a, /} - the color, spin and 
flavor indices. 

We use a minimal set of physical parameters, namely, the coupling constant, the 
scale of confinement and the quark masses: {a s = g 2 / 'in , A, m u d, m s , m c , nib} ■ Hereby, 
we do not distinct the masses of u and s quarks and neglect the superheavy t quark 
due to the absence of observed (tt) bound states. 

The matrix elements of hadron processes at large distance are in fact integrated 
characteristics of the quark and gluon propagators and the interaction vertices. Taking 
into account the correct global symmetry properties and their breaking by introducing 
additional physical parameters may be more important than the working out in detail 
(e.g., [9]). We admit that the bound states of quarks and gluons may be found as the 
solution of the Bethe-Salpeter equation in a variational-integral form [7]. One may 
guess that the solution is not too sensible on the details of integrands. 

We consider the effective quark and gluon propagators as follows: 

satf)= ^{tf + mll + T^M]}^,^ 

A 2 ( 2 A 2 ^ 

Df u {x) = 5 a %„ j^y 2 exp |- ^ j = 5 a \ u D(x) , (2) 

where p = p^^, m - the quark mass and < w(m) < 1. These entire analytic functions 
in Euclidean metric represent more realistic and reasonable extensions to model prop- 
agators used in our earlier investigations [5, 6, 7] and allow us to estimate the meson 
and glueball spectra with sufficient accuracy by avoiding cumbersome calculations. 
Consider the partition function 

z = JJs^s^ y^expj-^s- 1 ^) -i(0/rV) - g^r^cf)) - #(#f)J . (3) 
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We suppose, that the coupling a s is relatively weak (in Section 3.1 we show that the 
coupling is indeed small). Then, we can restrict the consideration within the one-gluon 
exchange approximation and take the Gaussian path integration over 0-variable: 

y^exp j-I(0irV) - g(m<f>) - g(<M>F)} = exp {-£ 2 [*, M>] - C G ) , (4) 
where, the terms corresponding to the two-quark and two-gluon bound states read 

C 2 = y Js<j> e- Cs ^ ((*rtt)£>(tfr*)) , C G = J6(f> e" £sM . (5) 

3 Mesons 

Describe shortly the important steps of our approach on the example of the two-quark 
bound states, more details can be found in [5, 6]. First, we allocate the one-gluon 
exchange between quark currents, go to the relative co-ordinates in the " center-of-mass" 
system and introduce the relative mass of quark \ij — mf/(mi + m 2 ). The latter is 
important for mesons consisting of different quark flavors. The next step is to perform 
a Fierz transformation to get the colorless bilocal quark currents. Then, introduce an 
normal basis {U(x)}. The appropriate diagonalization of £ 2 on the colorless quark 
currents should be done by using {U(x)}. Use a Gaussian path-integral representation 
by introducing auxiliary meson fields Bjj and apply the 'Hadronization Ansatz' to 
identify Bj^(x) with meson fields carrying the quantum numbers M = {«///'}, where 
J - the spin and / - the quark flavor. 

Below we consider the pseudoscalar (P : J PC = CT + ) and vector (V : J PC = 1 ) 
mesons, the most established sectors of hadron spectroscopy. 

The partition function for mesons reads: 

Z N = JDB X exp {~\{B N [1 + <? 2 Tr(W S V N S)\ B N ) + W^B*]} , (6) 



where VV = ^jjdy U(y)y D(y) exp{y/2f d} is a vertice function, J = {S, P,V, A,T} 
and Tj = {1,1^,1^^,1^^,(7^}. We use a Euclidean metric, with: {7 M ,7„} = 25^; 
lit = 7/ii a b = Z)i=i °A- For a timelike vector p^, p 2 < 0. 

The final-state interaction between mesons is described by W/[.B/v] = 0[B^\. 

The diagonalization of the quadratic form on the orthonormal system {Ujy} is 
nothing else but the solution of the Bethe-Salpeter equation (in the one-gluon approx- 
imation): 

1 + g 2 Tr(y M SV M ,S) = JJ dxdy U{x) jl + g 2 ^D{x) J e~ ik ^ 



TV 



(2tt)< 

TjS (k + ^p) Tj,S (k - /xap)] U(y) = 5 NN '[l + \ M (-p 2 )] ■ (7) 



Then, the meson mass is derived from equation: 

1 + A^(M^) = 0. (8) 
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3.1 Ground States 

We should solve the Bethe-Salpeter equation (7) with sufficient accuracy for the meson 
ground states. The polarization kernel is real and symmetric that allows us to find a 
simple variational solution to this problem. 

According to (7) it is reasonable to choose a normalized trial function [5, 6, 7]: 

Ui{^}{x, a) ~ T t{ll} (x)yjD(x) e" aA ^ , f dx |^m( x ' a )f = 1 > ( 9 ) 

M 

where i is the orbital quantum number, (x) is the four- dimensional spherical har- 
monic and a > is a parameter. 

Substituting (9) into (7) we obtain a variational equation for the meson mass 

= 1 + m&xJ2jjdxdy U i{fl} (x,a) g 2 U p (x,y) U i{fl} (y,a) , p 2 = -M} . (10) 
M 
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Figure 1: Evolution of the pseudoscalar meson mass with quark mass m u . 

Below we consider £ = 0, the pseudoscalar (J PC = 0~ + ) and vector (1 ) mesons, 
the most established sectors of hadron spectroscopy. 
For the ground state mass M = Mo we obtain: 
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a. s A 4 



M 2 {p\ + ill) - ml - ml 



487rmfm2 

bM 2 ^-^) 2 



2A 2 



max 

1/3<6<1 



(56-l)(l-6) 



1 2 



• exp 



4A 2 



M 2 ( , , , , (to - Ai 2 ) 2 ~ 

2b pj + — I K & + b(2- b P] ) ^ 1 



mim 2 



(11) 



where u(z) = 1/(1 + z 2 /4) and 

PJ = {1,1/2}, XJ = {+1,-1} for J = {P,V}. 

Equation (11) is nontrivial, the behaviour of the quark propagator S ~ 1/m 2 
imposes a restriction on the value of the quark mass from below. Particularly, the 
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evolution of the pseudoscalar meson mass Mp(m u , m u ) with "constituent-quark" mass 
m u at fixed values of {a s ,A} is depicted in Fig.l. The optimal value of m u is fixed 
from the relation Mp(m u ,m u ) = M n = 138 MeV. 

By comparing our estimates with the observed meson masses [8] we have found the 
following optimal values of free parameters [7]: 

a s = 0.186 , A = 730 MeV , m ud = 170 MeV , 
m s = 188 MeV , m c = 646 MeV , m b = 4221 MeV . (12) 
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Table 1: Estimated ground-state masses Mp and My (in MeV). 

Our estimates for the pseudoscalar and vector meson masses in the ground state 
(Tab.l and Fig.2) compared with experimental data show that the relative error is 
small (~ 1 -T- 5 per cent) in a wide range of mass ~ 140 -j- 9500 MeV. 
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Figure 2: The estimated meson masses (dots) and the experimental data (dashes). 
We note: 

1. The set of optimal parameters (12) is not unique, another set also gives a reason- 
able ground-state spectra for P and V mesons. But we choose these values by 
taking into account other important hadron data such as the orbital and radial 
excitations of mesons as well as the lowest-state glueball mass. 



5 



2. The obtained value a s = 0.186 is in agreement with the latest experimental data 
for the QCD coupling value a® ~ 0.10 -j- 0.35 on the hadron distance [8] and 
with the prediction of the quenched theory af 1 ^ 11 ~ 0.195 [10]. Moreover, this 
relatively weak coupling value justifies the use of the one-gluon exchange mode 
in our consideration. A possible use of a "sliding" interaction strength a s (M) 
depending on the mass scale (see, e.g., [11, 12, 10]) will be discussed later. 

3. The obtained quark masses fit an hierarchy: m u ^ < m s < m c < m^. 

4. Our model provides the S'?7(3)-symmetry breaking: M n ^ Mr. 

5. The to and $ mesons considered as mixed states of uu + dd and ss pairs with a 
priori unknown mixing angle By. The optimal value is found: 9y = 9°. 

6. The [/(l)-splitting is explained, i.e., we obtain a large mass difference M n /M v « 
0.25 between pseudoscalar isovector and isoscalar mesons, while M p jM^ « 1 for 
the vector isovector and isoscalar. 

7. The dependence of the meson masses M(m\+m2) on the sum of two " constituent" 
quark masses (Fig. 3) is nontrivial. 

8. A rough estimate of the hadronization region 

2 Jd 4 x x 2 D(x) 8 / 1 \ 2 
Thadr ~ J¥x x D{x) ~ A 2 ~ U50MeV/ 

shows that the confinement scale r con f ~ 1/A is comparable with rhadr- 
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Figure 3: The estimated masses for the P and V mesons versa the sum of quark masses. 
The dashed straight lines depict the plain sum mi + mi- 

Below, we extend our consideration to the orbital (£ > 0) and radial (n r > 0) 
excitations of mesons. 
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4 Meson Regge Trajectories with n r = 



The orbital excitations take place in larger distances and should be less sensitive to 
the short-range details of the propagators. Therefore, correct description of the meson 
Regge trajectories can serve as an additional test ground for our model. 

Below we concentrate mostly on the orbital excitations with zero radial number 
n r = for which experimental data exist for all states with £ < 5. 

Particularly, by substituting (9) into (7) we obtain the mass equation: 



VTLr, 



where u(z) 



ASumlml ¥(£ + 1) 



exp 



2A 2 



Ml 



1 



max 

1/3<c<1 



2c Pj + 



(5c-l)(l-c) 



1 2 



d 



° 2 Tc r exp 



(mi + m 2 ) 2 

cMf{m\ — m 2 ) 2 
4A 2 (mi + m 2 ) 2 



Ml 



A 2 (mi + m 2) 

fmi\ (m 2 



' c(2-c Pj ) , 
mi m 2 H — mi 



m 2 



mi 7^2 

-- 1/(1 + z 2 /A) and pj = {1,1/2}, Xj = {+1, -1} for J = {P,V} 



(13) 
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Figure 4: The estimated mass Mj (dots) for the p- and i^*-meson excitations. 

Substituting the optimal parameters (12) into (13) we have estimated the masses of 
excitations for the p-meson and if-meson families. Our estimates (dots) plotted versa 
the quantum number J = £ + s (s = {0, 1} - the spin of the (qq 1 ) pair) is given in 
Fig. 4 in comparison with the experimental data (boxes) [8]. 

As is expected, our model describes acceptably the J-trajectories of the p-meson 
and iT-meson excitations. 

For large £ 3> 1 the asymptotical behaviour of Mf reads 

Mf ~ 4(mi + m 2 f + 4A 2 (£ + 2) In ( 2 = M 2 + £ ■ Vl e . (14) 
We note that for large J ^> 1: 

i) The Regge trajectories of K* and p mesons are linear and similar. 

ii) The intercept Mq is positive and depends on quark masses and A. 
ii) The slope fit sensitively depends only on A. 
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5 Meson Radial Excitations n r > 



A nice example of the meson radial excitations is the "charmonium" sector. Experi- 
mentally well established is the J/^-family (cc) with J PC = 1 . 

We expand the radial part of the normalized trial function as follows: 

U n (x) ~ L^ ,2) (2bx 2 ) exp [-A 2 (bx 2 + cxl)) , J dx \U n (x)\ 2 = 1 , x = {x ,x} , 

where {b, c} - dimensionless positive parameters, n - the radial quantum number and 
(z) is the generalized Laguerre polynomial. 
Taking into account the orthogonality of L^l 2 ^ (z) on the interval (0, oo) with respect 



to the weight function w(z) 
spectrum: 



,1/2 



e z we arrive in the equation for the radial-excitation 
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d 2n 
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n 
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dt n ds n 



3uv 
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n\(n + l)\ a,b 



+ ^ + A -^(l-co 2 (m/A)) 



[(1 -t)(l - s) {uv + u + v)}*/ 2 



(15) 



t=0, s=0 



where 



u 



B + 



(2B-l)s 



B + 



{2B- l)t 
(l-t) 



<A<l<B<oo. 



The numerical solutions of Eq.(15) for M n (cc) are plotted in Fig. 5. We see that there 
exists an expressed convexity for n < 3 and for higher mass members the n-trajectory 
is approximately linear. By analysing (15) for large n we obtain the solution: 



M 2 (m, A) = M 2 (m, A) + n ■ Q n (A) . 



(16) 



For n ^> 1: 

i) The radial Regge trajectories are linear. 

ii) The intercept Mq depends on both m and A. 
ii) The slope Vt n depends only on A. 

Note, for the 1 radial excitations of charmonium, we obtain a correct (convex, 
asymptotically linear) behaviour, but underestimate the experimental data. For larger 
A ~ 800 + 850 MeV our estimate approaches the results in [8] . 



6 Glueball Lowest State 

In contrast to the meson sector, the experimental status of the glueball is not clear. 
The signatures naively expected for glueballs in the experiment are: the absence in any 
(qq) nonets, an enhanced production in gluon-rich channels of central productions and 
radiative decays, a decay branching fraction incompatible with two-quark states and 
the reduced couplings. 

There are predictions expecting non-gg scalar objects, like glueballs and multiquark 
states in the mass range ~ 1600 -=-1800 MeV [13, 14]. Recent lattice calculations, QCD 
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Figure 5: The spectrum of the J/^-family radial excitations for A = 730 MeV. 



sum rules, "tube" and constituent glue models agree that the lightest glueballs have 
quantum numbers J PC = ++ and 2 ++ . 

A QCD sum rule [15] and the K-matrix analysis [16] predict the scalar glueball 
mass Mq = 1-25 ± 0.2 GeV reducing down the quantized knot soliton model result 
[17] and the earlier quenched lattice estimates Mq — 1550 MeV [18, 19, 20]. However, 
most recent quenched lattice estimate with improved action favors a larger mass close 
to M G = 1710 ± 50 ± 58 MeV [21]. 

It is evident that the structure of QCD vacuum plays the main role in the formation 
of glueballs. We suppose that the lowest state of the glueball may be reasonably 
described by the " gluon-gluon" bound state J PC = ++ in our model. 

The Gaussian character of the gluon propagator (2) allows us explicitly solve the 
equation for the glueball mass and the result reads: 

M^2A^(^), avr = M^t, (17) 

Note, the glueball mass (17) depends on a s in a nonperturbative way and the 
physical restriction M G > leads to an upper bound: a s < a upp ~ 3.24124. 

Substituting the optimal parameters (a s = 0.186, A = 730 MeV) obtained for the 
meson sector into (17) we estimate the lowest-state scalar glueball mass 

M G = 1745 MeV . (18) 

Our estimate (18) does not contradict the latest predictions expecting glueballs in 
the mass range ~ 1600 1800 MeV [13, 14]. 

In conclusion, we have considered a relativistic quantum field model of interact- 
ing quarks and gluons under the analytic confinement. We have used only physical 
parameters (the quark masses, the coupling constant and the confinement scale) for 
the model and solved the Bethe-Salpeter equation in the ladder approximation for the 
hadron bound states. The use of one-gluon exchange mode is justified by the estimated 
small value a s = 0.186 of the coupling constant that is also in agreement with both the 
prediction of the quenched theory ot quenc h e d ~ 0.195 [10] and the latest experimental 
data for the QCD running coupling a s ~ 0.10 -j- 0.35. 
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Our approach does not require the "flux tube" -type confinement for the light and 
heavy mesons. 

Within a simple relativistic model with reasonable forms of the quark and gluon 
propagators we describe correctly: 

- the pseudoscalar and vector meson masses in the ground state and orbital excita- 
tions. The relative error is small in a wide range of mass 140 MeV 4- 9.46 GeV, 

- the quark mass hierarchy: m u ^ < m s < m c < mi,, 

- the SU(3)-symmetry breaking: Mk ^ M n , 

- the so-called "[/(l)-splitting": M n < M v while M p « M u , 

- the approximately linear radial and orbital Regge trajectories for the pseudoscalar 
and vector mesons, 

- the lowest state glueball mass close to ~ 1750 MeV. 

The author thanks V.V.Burov, G.V.Efimov, S.B.Gerasimov, E.Klempt andW.Oelert 
for useful discussions and comments. 
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